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ABSTRACT 



We study the collimation of relativistic magnetohydrodynamic jets by the pressure 
of an ambient medium, in the limit where the jet interior loses causal contact with its 
surroundings. This follows up a hydrodynamic study in a previous paper, adding the 
effects of a toroidal magnetic field threading the jet. As the ultrarelativistic jet en- 
counters an ambient medium with a pressure profile with a radial scaling of p oc r~ v 
where 2 < r\ < 4, it loses causal contact with its surroundings and forms a bound- 
ary layer with a large pressure gradient. By constructing self-similar solutions to the 
fluid equations within this boundary layer, we examine the structure of this layer as 
a function of the external pressure profile. We show that the boundary layer always 
becomes magnetically dominated far from the source, and that in the magnetic limit, 
physical self-similar solutions are admitted in which the total pressure within the layer 
decreases linearly with distance from the contact discontinuity inward. These solu- 
tions suggest a 'hollow cone' behavior of the jet, with the boundary layer thickness 
prescribed by the value of rj. In contrast to the hydrodynamical case, however, the 
boundary layer contains an asymptotically vanishing fraction of the jet energy flux. 
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1 INTRODUCTION 

The outflows from ac tive galactic nuc l e i (AG N) are thought to 
be highly r elativistic (Begelma n et al.lll.984) and highly colli- 
mated (e.g. ljorstad et aU 2005), but the cause of this collima- 
tion is uncertain. 

Because jet-launching is g enerally believed to be 
electromagnetically driven (e .g. [Blandford & Pavnel Il982t 
Contopoul os & Lovelace! [T9941) . one of the most commonly- 
accepted explanations for the observed collimation is that 
jets are threaded with magnetic fields that cau s e collima- 
tion via magnetic tension (e.g. Benfordl Il9"78fc iBegelmanl 
1995). Supporting this theory, it has been demonstrated 
that both relativistic and non-relativistic hydromagnetic out- 
flows must eventually becom e collimated dChiueh et aljl99 ll 
iHewa erts & Norman 1989]). For magnetic fields acting alone, 
howev er, collimation will only happen on extremely lar ge 
scales feichledll993tlBegelman & Lilll994lBegelmarJ[l995l) . 

To cause jets to collimate on reasonable scales, there 
must be an additional mechanism at work. A logical culprit 
is confinement by the pressure of an external medium. Pres- 
sure confinement has been demonstrated to act effectively on 



its own (e.g. Levinson & E ichler 20001 iBromberg & Levinsonl 
l2007l:lKohler et alj2012ir and accretion disk winds surround- 
ing an AGN provide an ideal ambient medium to help to col- 
limate the jet. 

There have been many numerical studies of magne- 
tized jets (e.g. iKomissarovi I1999L lHawlev & Krolikl 2006; 
iBeckwith etal]|2008l: iMcKinnev & Blandfordll2009l) . with the 
goal of forming a self-consistent description of the jet- 
launching and collimation mechanisms. These numerical sim- 
ulations have several restrictions, however, one of which be- 
ing that the boundary of the jet, rather than having its shape 
determine d by pressure balance, is genera l ly treated as a rigid 
wall (e.g. iKomissarov et all 120071 120091: iKomissarovl l201lh 
iTchekhovskov et alj I2O1O0 . This construction doesn't allow 
the ambient pressure to affect collimation of the jet. 

Treatments that do include effects of the external 
medium commonly focus on describing jets that remain i n 
causal contact (e.g. [Zakamska et alj|2008t ILvubarskvlbOllI) . 
As an ultrarelativistic jet expands into an ambient medium 
with a pressure profile p oc r^ 71 , it will eventually lose causal 
contact if 77 > 2 and the opening angle is greater than 1 /V, 
where T is the bulk Lorentz factor of the fluid. Observations 
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of gamma-ray bursts indicate that these relati vistic j ets largely 
have opening angles grea ter than 1/r (e.g. |Piranll2004l: see 
Tchekhovskoy et al. 2010 for discussion), and AGN outflows 
with large Lorentz factors may similarly be causally discon- 
nected; thus the poorly-studied regime of a jet that has lost 
causal contact is of physi cal interest. 

In a previous paper dKohler et al.l2012i hereafter KBB 12) 
we developed a model describing the recollimation boundary 
layer of a purely hydrodynamic, "hot" (pressure-dominated) 
jet with an ultrarelativistic equation of state. In this model, 
we assumed that the pressure outside the jet decreases with 
r so rapidly that the jet interior loses causal contact with 
its boundary, resulting in a shocked boundary layer form- 
ing within the jet. Though the jet interior is causally discon- 
nected, the boundary layer is nevertheless narrow enough to 
remain in causal contact itself. Assuming self-similarity as a 
function of r, we calculated how the transverse structure of 
the jet boundary layer depends on the value of 77 in the exter- 
nal pressure profile. 

We now expand this work to include, in addition to col- 
limation by the external medium, the effects of a magnetic 
field within the jet. We include only a toroidal field, as it is the 
toroidal field that domin ates the dynamics at large radii, far 
outside the light cyl inder ( Begelma n et aLll984t lContopoulos 
1995; Beskin 2009). In Section 2, we first demonstrate that 
seeding a jet with a magnetic field at the base will always 
cause it to become magnetically dominated at large radii. We 
then find self-similar solutions for the boundary layer of the 
jet in the limit of magnetic dominance. In Section 3 we discuss 
the results, and in Section 4 we conclude. 



2 SELF-SIMILAR TREATMENT OF THE MAGNETIZED 
BOUNDARY LAYER 

We use spherical coordinates to model a hot, ultrarelativistic 
jet that is symmetric about the z-axis and has approximately 
radial streamlines - an approximation justified because the 
jet interior is causally disconnected from the environment. We 
assume that the jet is injected from a point source with steady 
flow, and we examine the jet in its steady-state configuration. 

We focus on modeling the boundary layer of jet material 
that forms at the interface between the jet and the station- 
ary ambient medium. This layer is bounded on the inside by 
a shock front or a rarefaction front, and on the outside by 
a contact discontinuity. There is no mass flux across the con- 
tact discontinuity, and the pressure must be matched on either 
side of it. The physics of the ambient medium is wrapped into 
the external pressure profile p e . 

We adopt a pressure profile for the ambient medium of 
p e oc r~'', fixing the pressure external to the jet to be depen- 
dent only on the parameter 77. We focus on the case of ambient 
pressure where 2 < r? < 4, as in our work in KBB12 or in, e.g., 
iBromberg & Levinsonl d2007l) . because the jet interior is out of 
causal contact with the exterior for this range in r\. Physically, 
this pressure profile range could describe a confining wind, 
an accretion flow, a disk coron a, or even a stellar e nvelope in 
a GRB collapsar model (e.g. lBromberg etailboill) . 

As in KBB 12, we assume that the opening angle of the jet 
is much greater than 1 /V, such that causal contact has been 
lost. We construct a boundary layer that remains in causal 
contact, such that its thickness is of order A6 ~ 1/r. The 



boundary layer is thus very thin compared to the width of 
the jet. The radius of curvature of the jet is then much larger 
than the width of the boundary layer, allowing us to treat the 
curvature as a small effect. 

We begin with the equation s for relat i vistic MHD 
(RMHD) in flat spacetime ( e.g. iDixonl Il978t iKomissarovl 
Il999fc IZakamska et alj [2008), including the effects of a 
toroidal magnetic field within the jet. We ignore rotation since 
our regime of interest is far outside the light cylinder and, 
indeed, far outside the fast magnetosonic surface, rendering 
rotation effects unimportant. We again assume an ultrarel- 
ativistic equation of state such that the total proper energy 
density is given by e — p + 3p « 3p, where p and p are, re- 
spectively, the proper rest mass density and the pressure of 
the fluid within the boundary layer. 

The continuity equation remains unchanged with the ad- 
dition of a magnetic field, 

v ■ ( P pr) = 0, (i) 

where /3 and T are the velocity (/3 = v/c) and the bulk 
Lorentz factor of the fluid within the boundary layer. Assum- 
ing an energy equation of p oc p 4,/3 and taking 6 ~ constant, 
the continuity equation becomes 

^9 y /4 id m 

r 2 or r 06 

in spherical coordinates. 

Denoting the observer-frame magnetic field within the jet 
boundary layer as B, and using the ideal MHD condition to 
express the observer-frame electric field as E = — v x B, the 
momentum equation can now be written as 

(4pT 2 + B 2 )((3 ■ V)/3 + V(p + \b 2 Y- 2 ) 

-s[v-(Br- 2 )]-r- 2 (s-v)B = o, (3) 

where the parallel and perpendicular components of the mo- 
mentum equation are obtained by taking the respective vector 
dot and cross product of f3 with Eq ^ . 

Finally, we add the equation for flux freezing, 

V x (3 x B) = 0. (4) 

Writing this in spherical coordinates, and assuming a toroidal 
magnetic field B = B(r, 6)<p, we have 

§- r {rp r B) + ^- e {p B B) = Q. (5) 

The full form of these four equations admit self-similar 
solutions only in the case where r\ = A, which corresponds 
to a quasi-monopole flow with no collimation and therefore 
no distinct boundary layer. For all other values of r\, the full 
equations allow only for trivial solutions due to overconstraint 
of the system. We now demonstrate, however, that it is not 
appropriate to use these equations in their full form. They 
should instead be examined in the magnetically-dominated 
limit - where they do admit non-trivial self-similar solutions. 

2.1 Demonstration of Asymptotic Magnetic Dominance 

Suppose that the magnetic field takes the form 
B = rp 3/4 Tf(r,6) where / is some function of r and 
9. Inserting this into Eq (0, we obtain 

V ■ (p 3/4 r/3,f) = (6) 
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which can be expanded and combined with Eq ([XJ) to show 
that 



(7) 



indicating that the function / must be constant along stream- 
lines. 

Because the contact discontinuity is a streamline, we can 
therefore s tate that B oc rp 3 ^Y a long the contact discontinu- 
ity. As in lZakamska et all fe008l) . we now define a magneti- 
zation parameter Pb as the ratio of magnetic to gas pressure 
(the inverse of the usual plasma beta) . Along the contact dis- 
continuity this parameter is thus given by 



a & 2 1/2 

/3b = — ^ <xr p ' . 



pY 2 



(8) 



Because pressure must be matched across the contact dis- 
continuity, the external pressure p e oc r~ v must be balanced 
at that point by the total internal pressure within the bound- 
ary layer, ptot = p + B 2 /Y 2 . We now examine /3b in two 
extreme cases: the limit where the pressure balance at the 
contact discontinuity is supplied solely by the gas pressure 
within the boundary layer, and the limit where the balance is 
supplied solely by the magnetic pressure within the layer. 

In the gas pressure-dominated case, the internal gas pres- 
sure p is equivalent to the external pressure p e , and must 
therefore scale in the same way, such that p oc r~ v . Applying 
this to Eq $8^ demonstrates that in this case, /3b oc r^ 71 ^ 2 . 
Thus, for 2 < r\ < 4, /3b scales as r to some positive power. 

In the magnetically-dominated case, the external pres- 
sure is balanced by the magnetic pressure such that B 2 /Y 2 oc 
r~ v . This scaling implies that the internal gas pressure is 
given by p oc r -W 3 M+ 2 \ and we obtain /3 B oc r (4 "' ,)/3 . 
Again, for 2 < 77 < 4, /3b scales as r to some positive power. 

Thus we see that in both extreme cases, [3 B grows with 
increasing r. This suggests that no matter how small a mag- 
netic field the jet is seeded with, the boundary layer will 
eventually become magnetically-dominated far from the jet 
source. With this in mind, we now repeat the calculations 
performed in KBB12 with the inclusion of a toroidal magnetic 
field, specifically in the limit where {3b 2> 1. 



2.2 Solutions in the Magnetic-Dominance Limit 

We first rederive the fluid equations in Section 2.1 in the limit 
where (3b S> 1. Continuity and flux freezing are unchanged, 
but terms in the momentum equation containing I //3b are 
negligible in this limit. 

We now make scaling arguments as in KBB12: we as- 
sume Pe is of order 1 /Y, since this is the maximum transverse 
speed that can be achieved without a shock forming, and f3 r 
is of order one. With this characteristic scale, 



Ex- 



dp a dp B / dB a dB\ 



(9) 

02 ( dp e d/3 e 1 8T\ dp B dB 

(10) 

We now attempt to construct self-similar solutions in the fol- 
lowing fashion: 



P = 9(0r~ 
B = b(£)r- 



-77/2 



Pe = H0r- x , 
p = a(£)r~ a , 



(11) 



such that the external gas pressure is matched by the internal 
magnetic pressure, B 2 /Y 2 oc r~ v . In these solutions x and 
a are constant free parameters describing the radial scaling, 
and g, h, b and a are functions of a similarity variable £ (as in 
KBB12) that describes the distance from the contact disconti- 
nuity, normalized by the expected scale of the boundary layer, 
£ oc (0 C — 0) /AO. The angular thickness of the boundary layer 
is expected to scale as A9 = 1/T C , such that £ oc r x (8 c — 8), 
where 8 C = 8 c (r) is the location of the contact discontinuity. 

The fact that the streamlines at 8 C must be parallel to the 
contact discontinuity, requiring that Pg(8 c ) = rd8 c /dr, yields 
the further constraint that 



dr 



= hor 



-(*+!) 



(12) 



where ho = h(£ = 0) is a negative constant for collimat- 
ing solutions. While the flow is very nearly radial, this ex- 
pression describes the small deviation of the flow lines re- 
sulting from subtle collimation. Choosing the proportionality 
constant such that £ is defined as 

Z = ~r*(6 c -9) (13) 
no 

absorbs the boundary condition into the similarity variable 
and ensures collimating solutions (such that ho < 0) . 

We now recast the fluid equations in terms of these func- 
tions. The continuity and flux-freezing equations are fully self- 
similar and become, respectively, 

_ , h \ ( . „ „ ,h'' 
x£ - 1 + 



^ 9 



ho 



x£ - 1 + 



- ( 4a; + 8 - 3a + 4— J = 
(14) 



ho 



+ x 



2 +1+ h- 



= 0, 
(15) 



where primes denote differentiation with respect to £. 

Now examining the parallel component of the momen- 
tum equation, one can see that the radial scaling does not 
automatically vanish: 

b 2 9 2 



l+x -i + - b 



ho 



+r" 



ho 



(16) 



pressing p r in terms of pe and Y and employing the fact that Because we are specifically examining the regime where Pb oc 



Pe + Y~ 2 < 1, we have p r « 1 - §($ + Y~ 2 ). Using these B 2 /(pY 2 



3> 1, however, we can assume that r v a 



scalings and keeping terms only to lowest order, the parallel 
and perpendicular components of the momentum equation 



will be very small for large r, rendering the second term neg- 
ligible. 
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Using this logic, the parallel and perpendicular compo- 
nents of the momentum equation become, respectively, 



a 



V b' ( , , h 



h(l-x)+h'[xt-l+£ +f ^- + T 



b' 



= 



0. 



(17) 
(18) 



Assuming that b, g, a and h are finite and non-zero, Eqs 
(1 1 51) and 1 11 71 1 imply that h' = 0. This assumption results in 
Eq 015D becoming 



6' _ x - 77/2 + 1 

7T ~~ ~x~i 



(19) 



and yields the following general solutions for the functions 
describing the transverse behavior within the boundary layer: 



ho 

££-(2+2x-r,)/2x 



a 
a 



± 
3a 



{2+2x-r))/x _ 2h-Qx(x — 1) 
2 + x - 11 



-,1/2 



2r, 



3.T 



(20) 
(21) 

(22) 
(23) 



where ^4 and J5 are constants of integration that are defined 
by the boundary conditions. 

To produce physical solutions, we examine the special 
case where we prevent b' from having a singularity at £ = 
by setting x — rj/2 + 1 = 0, implying that 6 is a constant. In 
terms of boundary conditions h ,g ,b , and a , which serve 
as scaling factors and allow us to determine our functions self- 
consistently, the physical solutions within the boundary layer 
are therefore 



ho 
b 



9 = -ho 



a = ao 



1/2 



4/3 



-(4-f?)£ 



2/3 



(24) 
(25) 

(26) 
(27) 



with the constraints x 



p = 3?~ 



1-17/2 



- 1 and a = + 2), such that 

„ _ „„-(2/3)0l+2) 



(28) 



A check for self-consistency shows that the density in the 
lab frame, given by p 3/ ^ 4 F, has no dependence upon 6 and 
scales as r~ 2 , as is expected for nearly radial flow. 

It should be noted that for the toroidal magnetic field and 
corresponding electric field to exist within the boundary layer, 
there must be a current distribution and charge distribution 
within the layer, and a current sheet and surface charge at the 
outer boundary where the magnetic and electric fields termi- 
nate. Calculating the current distribution within the jet from 
the solutions in Eq 1 128D , one can see that longitudinal current 
within the jet is conserved in the case of approximately radial 
streamlines, providing another check of self-consistency. 




Figure 1. Special-case solution for in terms of the pressure- 
profile parameter jj and boundary conditions go and ho . 



3 DISCUSSION OF RESULTS 

The first important result of this solution is that, in the ob- 
server's frame, both the magnetic field B and the transverse 
jet velocity ft have only radial dependence; they are constant 
across the boundary layer. This is in direct opposition to the 
results from the strictly hydrodynamic limit (see KBB12, Sec- 
tion 3.1), where ft decreases monotonically from the outside 
of the boundary layer inward for all cases where 2 < 77 < 4. 

Another significant point is that the solution for the mag- 
netic field has no dependence in either dimension on the 
parameter 77, meaning that the magnetic field that develops 
within the boundary layer is not affected by the pressure pro- 
file of the medium that the jet passes through. 

The Lorentz factor V, on the other hand, does have a 
radial dependence on r\. As expected, we see that the jet is 
accelerated as it propagates outward: V oc r'" 2_1 scales as a 
positive power of r. 

We now examine the pressure profile within the bound- 
ary layer as prescribed by this solution. Since we have demon- 
strated that the magnetically-dominated regime is the rele- 
vant regime in this problem, the total pressure within the 
layer is approximately described by the magnetic pressure 



B 2 , > , > 



£l -<«-»« 



(29) 



In KBB12 we demonstrated that in the strictly hydrody- 
namic case, the pressure monotonically decreases for all 77, 
but decreases linearly with £ only for the case where 77 = 8/3 
(KBB12, Section 3.1). In this specific case, the pressure drops 
to zero within the boundary layer, indicating that all of the jet 
material is piled in a thin boundary layer in a 'hollow cone' 
structure. As the value of 77 gets further from 8/3 in either 
direction, the pressure profile becomes less steep, implying 
that the boundary layer broadens and the structure of the jet 
becomes less like a hollow cone. 

The magnetic case is a little more difficult to interpret 
due to the unspecified boundary conditions go and ho in Eq 
I I29D , but we can make some qualitative observations. First, 
it is clear that the magnetic case parallels the hydrodynamic 
case in that all solutions for pressure (which, since magneti- 
cally dominated, scales as 1 /T 2 in the boundary layer) mono- 
tonically decrease with £ for 2 < 77 < 4. This means that the 
pressure is greatest at the contact discontinuity where it is 
matched with the external pressure, and it decreases inwards 
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across the layer (see Fig. [TJ, as is expected due to the colli- 
mation. 

A major contrast between the magnetic and hydrody- 
namic cases is that in the magnetic case it is true for all values 
of j) that the pressure decreases linearly and vanishes at a crit- 
ical value of f, in this case £ cr = {go/ho) 2 /(4 — rj). Thus any 
jet can form a hollow cone structure in the magnetic case, 
with a boundary layer of thickness A£ = £ cr . The impact of 
?7, then, is instead to determine the thickness of the boundary 
layer: jets injected into an ambient medium with a pressure 
profile parametrized by a low r) (a gradual decrease in ex- 
ternal pressure) will therefore develop the most pronounced 
hollow-cone structure, whereas jets injected into a medium 
that has a higher 77 (a steep decrease in external pressure) 
will exhibit a thicker boundary layer. 

The distinction between the magnetic and hydrodynamic 
cases is most evident, however, when examining the global 
energy constraints for these solutions, as in KBB12. In this 
case, we assume that all the power in the boundary layer is in 
the form of Poynting flux, given by L oc B 2 A, where A is the 
cross-sectional area of the boundary layer. As we are assuming 
a boundary layer of width A8 = 1 /F, the cross-sectional area 
of that ring is A oc r 2 /r, implying that the power within the 
boundary layer scales as L oc r 1 "''/ 2 . Thus for 2 < 77 < 4, 
the boundary layer contains less and less energy as one goes 
further out in r. 

This loss of energy is indicative of a unique situation 
in the magnetic case: rather than being bounded on the in- 
side by a shock front through which material is added, as in 
the hydrodynamic case, at large r the boundary layer instead 
appears to be bounded on the inside by a rarefaction front, 
through which material is leaving the layer and returning to 
the main jet. 

We can verify the position of this rarefaction front by cal- 
culating the location of the surface at which the fluid motion 
normal to the surface is equivalent to the local speed of sound. 
This condition can be stated as 



01 



■ft 



(30) 



where /3^ r is the velocity normal to the surface of the rarefac- 
tion front f r , and /3; is the local sound speed, which is given 
by the relativistic expression for the fast magnetosonic speed 
in the hot-jet limit, 



PI 



Ip + b 2 



(31) 



4p + B 2 ' 

This calculation yields an approximate position for the 
rarefaction front of 



fa -2) 



3/2 



2ao 



3/4 



|0 ^ /2 r (,-4)/4 ] (32) 

ho , 



- ^ (4 - ,7)5/2 

indicating that the rarefaction front occurs at a location just 
before the pressure within the boundary layer drops to zero, 
allowing matching across the front to the conditions in the 
interior of the jet. At large r, the position of the rarefaction 
front asymptotes to the location of the pressure zero-point. 

Looking at how the Lorentz factor T scales with r along 
the rarefaction front, we see that T oc r 3v/6 - 1/2 . We know 
that the Lorentz factor scales as T oc r'"' 2 " 1 at the contact 
discontinuity, and we expect it to scale as T oc r i n the jet in- 
terior to be consistent with free expansion (see also Lyubarsky 



120111) . Thus, for 2 < rj < 4, the scaling of r at the rarefaction 
front is consistent with an intermediate acceleration between 
the two, providing a smooth transition between the boundary 
layer and the jet interior. 

More insight into the boundary-layer energy loss can be 
attained by examining the position of the contact discontinu- 
ity {8 C ) and the inner boundary (8 S = 8{g = 0)) as a function 
of radius. Using the solutions found in Eq H28D . the width of 
the boundary layer is given by 



9o /ho i- 

"~J r 

4 — rj 



1/3 



(33) 



indicating that the width is decreasing in r. 

Furthermore, integrating Eq ( 1121 ) using the above equa- 
tion allows us to find the form of the inner and outer bound- 
aries individually: 



9o 



ho 



r,/2 - 1 
ho 



n 



-07/2-1) 



4 - 77 

*?/2-l 



-("7/2-1) 



(34) 



(35) 



From this, we can see that the outer boundary is collimating, 
but less and less so with increasing r. The inner boundary, on 
the other hand, is decollimating when the condition 



77/2- 



(36) 



is met, with the boundary conditions prescribed by match- 
ing across the rarefaction front to solutions for the jet inte- 
rior (such as those described by Lvubarsky 2011). When this 
condition is satisfied, the boundary layer intercepts fewer and 
fewer streamlines of new material from within the jet interior. 
The decollimation, too, weakens with increasing r, with the 
inner and outer boundaries meeting asymptotically. 

These results suggest that the boundary layer of a mag- 
netically dominated jet decreases in width with increasing dis- 
tance from the source, and contains a decreasing amount of 
energy as material leaves the boundary layer across the rar- 
efaction front to rejoin the jet interior. Nonetheless, a sharp 
pressure gradient is maintained across the layer, insulating 
the interior of the jet from the external medium. 



4 CONCLUSION 

We have evaluated the structure of a boundary layer within a 
hot, magnetohydrodynamic jet with an ultrarelativistic equa- 
tion of state. We assumed that the jet as a whole is causally 
disconnected from its surroundings, but the boundary layer 
is thin and therefore in causal contact. We examined the im- 
pact on jet collimation of a toroidal magnetic field within the 
jet as well as an ambient medium with a pressure profile of 
p oc r - '', with 2 < 77 < 4. 

We first demonstrated that the basic RMHD equations 
can be used to show that any jet boundary layer seeded with 
a toroidal magnetic field at its base will eventually become 
magnetically dominated at large radii. 

We then constructed self-similar solutions for the bound- 
ary layer in the limit where the jet pressure is dominated by 
magnetic pressure. We found a special case of physical so- 
lutions where the jet pressure decreases linearly across the 
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boundary layer, dropping to zero at a location set by the 
boundary conditions and the value of the pressure profile pa- 
rameter rj. The boundary layer thickness is dependent upon 
the value of 77, with increasing 77 producing an increasingly 
wide boundary layer. 

We further found that the thickness of the boundary layer 
decreases with radial distance, and the boundary layer con- 
tains a decreasing amount of the jet energy. This suggests that 
the addition of a magnetic field fundamentally changes the jet 
at large radius: whereas in the hydrodynamic case the inner 
boundary of the layer is a shock front through which material 
enters the layer, in the magnetohydrodynamic case the layer 
is bounded on the inside by a rarefaction front through which 
material leaves the boundary layer and rejoins the interior of 
the jet. 

We found the position of this rarefaction front to occur 
just inside the layer from the location where the pressure van- 
ishes, providing a smooth transition between the boundary 
layer and the jet interior and allowing for matching across 
the rarefaction front to the conditions in the jet interior. This 
matching would prescribe the values of the boundary condi- 
tions within the layer, and could potentially yield a solution 
where the rarefaction front is gradually decollimating, inter- 
cepting fewer and fewer streamlines as radial distance from 
the source increases. 

In spite of the thinning of the boundary layer with ra- 
dial distance, a sharp pressure gradient is nonetheless main- 
tained across the layer, causing it to function as an insulating 
buffer between the jet interior and the ambient medium. Un- 
like the hydrodynamic case, the solutions for the structure of 
a magnetized jet do not have clear observational implications. 
Though the boundary layer contains a decreasing amount of 
energy as one looks further from the source, the layer might 
nonetheless have a high emissivity, which could be observa- 
tionally important if the flow within the boundary layer is 
pointed along our line of sight. 

The results presented in this paper provide the premise 
for a more complete treatment in numerical simulations of 
the effects of the ambient medium on collimation, both by 
demonstrating the behavior of the jet when the outer wall is 
allowed to change its shape, and by providing models that 
can be used to assess the effects of numerical resolution on 
simulation outcomes. 

Ultimately, these results provide a foundation for future 
work examining energy dissipation in magnetized jets and the 
associated radiative observational signatures. 
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